Dynamical formation and manipulation of Majorana fermions in driven quantum wires 
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Controlling the dynamics of Majorana fermions (MF) subject to time- varying driving fields is of 
fundamental importance for the practical realization of topological quantum computing. In this 
work we study how it is possible to dynamically generate and maintain the topological phase in one- 
dimensional superconducting nanowires after the temporal variation of the Hamiltonian parameters. 
Remarkably we show that for a sudden quench the system can never relax towards a state exhibiting 
fully developed MF, independently of the initial and final Hamiltonians. Only for sufficiently slow 
protocols the system behaves adiabatically, and the topological phase can be reached. Finally we 
address the crucial question of how "adiabatic" a protocol must be in order to manipulate the MF 
inside the topological phase without deteriorating their Majorana character. 
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Introduction. — The enormous potential of Majorana 
fermions (MF) for implementing decoherence-free quan- 
tum computation[l| is greatly stimulating their search 
in solid state systems. Several different platforms have 
been proposed for realizing MF, like fractional quan- 
tum Hall states at filling factor 5/2, vortex cores of 
p + ip superconductors [2J, surfaces of topological in- 
sulators coupled to s-wave superconductors 3] , non- 
ccntrosymmetric superconductors I4j , and ferromagnetic 
Josephson junctions[5j. However, despite such intense 
theoretical activity, no experimental evidence of the ex- 
istence of MF in the above systems has been provided so 
far. A promising alternative that should drastically sim- 
plify the realization and detection of MF consists in form- 
ing one-dimensional (ID) heterostructures with semicon- 
ductors and conventional superconductors 0,01 • Here the 
p-wave superconductivity is simulated by means of spin- 
orbit (SO) interaction and the system can be driven into 
the topological (T) phase (where the MF appear) by ap- 
plying a suitable magnetic field. In addition MF occur- 
ring in 2D and 3D networks of these nanowires are par- 
ticularly attracting, since they can be adiabatically ma- 
nipulated by using tunable gates or by reorienting of the 
magnetic field 7, 8]. Very recently signatures of MF pres- 
ence in InSb and In As nanowires have been reported 0- 
[l2T ] . Quantum wires hosting MF could also be engineered 
in trapped ultracold fermionic atoms, by employing op- 
tical Raman transitions to generate effective spin-orbit 
coupling and Zeeman fields, and using the proximity ef- 
fect with a bulk molecular Bose-Einstein condensate (l3$. 
In these systems one can tune the parameters with high 
temporal precision and efficiently control the amount of 
disorder. Thus they offer a unique possibility to study in 
a very clean way the nonequilibrium dynamics of MF fol- 
lowing time-dependent perturbations, and in particular 
their formation when the system undergoes the topolog- 
ical phase transition. 

In this Letter we study the real-time dynamics of a 
ID quantum wire of finite length in contact with a su- 



perconductor and in presence of SO coupling. The sys- 
tem is driven out of equilibrium by varying the external 
magnetic field by means of different protocols. The for- 
mation of MF is monitored by calculating the evolution 
of the recently proposed Majorana order parameter 
Remarkably we show that for sudden variations of the 
Hamiltonian the system relaxes towards a nonthermal 
state that does not exhibits fully developed MF. Only for 
sufficiently slow protocols the system undergoes the topo- 
logical phase transition, over a timescale which increases 
with increasing length of the wire. Finally we show that 
the manipulation of MF inside the T phase must obey 
precise temporal constraints in order to preserve the Ma- 
jorana character. This is a crucial issue for the practical 
implementation of topological quantum computing. 

Model and formalism. — The Hamiltonian of the cou- 
pled wire of length C is given by0, HH 

H = [ dx^(x)(-dl/2m - fj, - iad x o y + V z a z )V(x) 
Jo 

+ A^ t (x)^(x) + h.c.), & = y,\,^l), (1) 

where ui are the Pauli matrices, ip a {x)^' annihilates (cre- 
ates) an electron of mass m and spin a at position x in 
the wire, \i is the chemical potential, A is the strength 
of the proximity pairing field, and a and V z are the am- 
plitudes of SO coupling and Zeeman field respectively. 
This model displays a topologically trivial (TT) phase 
for Q = /i 2 + A 2 - V? > and a T phase for Q < 0, 
the phase transition occurring at Q = 0@, @]. The exis- 
tence of MF in the T phase can be probed by considering 
the Majorana polarization 14j . defined as the anomalous 
local density of states 

1 f°° 

P(x,u)=—Tm / e iw *V2*({^M),^M)}), 
n Jo V 

(2) 

where the average is taken over the Hamiltonian ground- 
state |\&o)- As discussed in Refs. 14. ltij a good order 
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parameter to detect the T phase can be built from P as 



nC/2 

$ = / dxP(x, 0) 
Jo 



dxP(x,0). 



(3) 



C/2 



Indeed in the TT phase it holds $ = 0, while in the 
T phase only the MF contribute to the order parameter 
yielding $ = 1. The nonanalytic behavior of $ at Q = 
indicates the occurrence of the phase transition, while 
values of $ close to 1 signal MF not fully developed. It 
is worth observing that for noninteracting fermions [as 
in the case of Hamiltonian in Eq. (JTJ] the polarization 
takes the simpler form[I3| 

P(x, w) = 2 £ <J(w - e n )uW*(x)v^ (x), (4) 



where n denotes the n-th eigenstate of the system with 
eigenvalue e n and wavefunction (ttj. , i?j n ' , , ) ex- 
pressed in the basis (■)/>}, ip^, ^j) V*t)- 

The pair of (real) zero-energy MF supported by the 
Hamiltonian H in the Tphase are located at the two op- 
posite edges of the wire[6|,ljj and decay exponentially into 
the bulk. The corresponding wavefunction must obey the 
constraint u a — v„ (that ensures the particle/ antiparticle 
equivalence) and can be found by solving the auxiliary 
problem 6] 



-d%/2m- (j, + V z -77A 



r]A - ad x 



-dl/2m 



ad x 

■n-v a 



u t (x) 



(5) 

where rj = ±1. The two choices of r\ provide the MF 
located at the right and left boundary of the wire re- 
spectively. In the following we assume r\ = 1, since the 
calculation with r\ = — 1 follows the same line of reason- 
ing. The MF wavefunction is then obtained by imposing 
the ansatz 



u t (x) 



_fu f 



(6) 



with g(x) — e~ x '^[17|. The characteristic equation for £ 
has in principle 4 complex solutions for a given I. How- 
ever, in the T phase (i.e. for Q < 0) and for any 77 = ±1 
only a single solution is real and ensures the normalizabil- 
ity of the wavefunction Q. In addition it is straightfor- 
ward to verify that the allowed £ for 7/ = ±1 do coincide, 
as dictated by symmetry constraints. The coefficients U a 
are easily obtained and read 



(-1/2 



'1 + 



V z - 



a/e+A 



1-1/2 



(7) 
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From the above solution we obtain that the Majorana 
order parameter $ = 2r]£J2 a ^ s 1 m the left half- wire 
and —1 in the right half- wire, as it should beflif. 

Real-time evolution. — We now study the real-time 
evolution of $ after the variation of the external magnetic 
field V z according to different protocols. The explicit cal- 
culations are performed within the tight-binding version 
of the Hamiltonian in Eq. (TTJ), which reads 



H 



AT 

£ 

i=l 

a 



CtC i+1 +h.c.-( f x-v)ClC i 



~(iC}a y C i+1 + h.c.) + V z C\a z d 
+ A( Clt o 4 + h.c), Cj = (cl v cl), 



(8) 



with M = C/ a and a the lattice spacing. The map- 
ping between the parameters of the continuum and lat- 
tice models is discussed in Ref. [l8j]. In the follow- 
ing we express energies in units of the hopping v and 
times in units of 1/v. If V z — > V z (t) the Hamiltonian 
becomes explicitly time-dependent H —t H(t) and the 
dynamics of the system is addressed by propagating the 
equilibrium ground-state l^o) appearing in Eq. © ac- 
cording to | v I / o(^)) = Texp[— i J Q dsH(s)]\^fo), where T 
is the time-ordering operator. The problem is numer- 
ically solved by discretizing the time and calculating 
the evolution of j^o) within a time-stepping procedure 
|*o(tj)> ~ eaq>H-ff(*j)**]|*o(*j-i)>, where t 3 =J^t^5t 
being a small time step and j a positive integer [19l [20|. 
We have considered ramp-like switching protocols of du- 
ration t bringing the magnetic field from the initial 
value V z at t = [with corresponding Hamiltonian 
H(0) = Hi] to the final value V z [with Hamiltonian 
H (t > t) = Hf] which is maintained constant for t > r, 

i.e. V z {t) = e{r-t)[V^ + {vP-V z {i] )}t/T+e{t-T)V z U \ 
The order parameter is then extracted by following Ref. 

Since energy is not conserved during the temporal 
evolution, we calculate $(t) by integrating over x the in- 
stantaneous polarization P(x, E op t(t)), where E opi (t) is 
the minimum average of H(t) over the evolved cigenstates 
of Hi, denoted by \<j>„ (t))[2lj|. At every time the states 
corresponding to the mininum energy are always two, 
with same value of |-E pt(i)|> but with opposite signs due 
to the symmetry of the problem. Therefore the two opti- 
mal \4>opt{t)) [with wavefunctions u^ pt \x , t) , v^ pt \x , t) 
to be inserted in Eq. (J5J] are the best approximation to 
the pair of MF that the system can provide at a given 
time. 

Quench dynamics. — If r = + we have a so-called 
sudden quench. We have studied two relevant situations, 
namely (i) the system initially in the TT phase and the 
value V z such that Q < for the final Hamiltonian 
(TT— >T quench), and (ii) Q < both for the initial and 
final systems (T — > T quench) . The first case serves to 
understand how the MF are dynamically formed, while 
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FIG. 1: $ as a function of time for a TT — > T quench with 



(0 



10~ b and V z 



(/) 



0.2 (left panel), and for a T -> T 



quench with Vr = 0.2 and V z u > = 0.3 (right panel). The 
rest of parameters are /i = 0, A = 0.1, ot = 0.3. The dashed 
line in the right panel represents $ggb calculated as in Eq. 
(fTTJll . $gge has been obtained with Af = 400, but its value 
depends very weakly on N '. 



the second case is crucial to investigate whether the MF 
maintain their character after a manipulation of the sys- 
tem parameters inside the T phase. In Fig. [T]we show the 
time evolution of the order parameter in the two cases, for 
different values of the size Af. Remarkably if the system 
is initially in the TT phase the MF do not form after the 
quench, and the order parameter remains close to zero 
for any Af, displaying temporal oscillations due to the 
finite size of the system (Fig. [T] left panel). If instead 
the system is initially in the T phase, we see that the 
MF of Hi are corrupted by the quench dynamics even 
if the condition Q < is maintained at all times, and 
$ approaches finite values significantly smaller than 1 
(Fig. [T] right panel). Also in this case $(f) displays oscil- 
lations which, however, tend to disappear by increasing 
Af. We have verified that similar results are also found 
by quenching other quantities like the strength of the 
pairing field, the SO coupling or the chemical potential 
(not show) . The above findings indicate that after a sud- 
den quench the initial ground-state |^o) does not relax 
to the ground-state of the quenched Hamiltonian 22 - 25| , 
and, more in important, that no single-particle eigen- 
state of Hi transforms into the MF of Hf. It has been 
conjectured 26] that some of the properties of the non- 
thermal state that develops after a sudden quench can 
be addressed within the so-called generalized Gibbs en- 
semble (GGE). The GGE permits to compute a class of 
long-time averages of of integrable quenched systems by 
means of the special density matrix |26j 



1 



Pgge 



J GGE 



- £„ A„x„ 



(9) 



where Zqge = Tr[e~ ^ n XnXn ] and I n are a set of in- 
tegrals of motion of the quenched Hamiltonian. The 
weights l n are fixed by the condition TrfpcGE^n] = 
(2n)t=o- Thus one can argue that the average of an 
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FIG. 2: <E>gge for a T — > T quench as a function of at 
fixed Vz — 0.2 and V z = 0.3 for both the continuum and 
the lattice models. The rest of parameters are the same as in 
Fig. [JJ We set a = 1 and hence m = l/i> 18]. For the lattice 
model we computed numerically Eq. ()10|l by diagonalizing 
H with Af = 200 (the result is however quite insensitive on 
the value of Af) while for the continuum model we used the 
analytic wavefunction in Eqs. ()6I7|) . 



observable O after the quench is obtained as Ogge = 
Tr[/9GGE0]- In the present case we choose X n as the 
eigenmode-occupations of the quenched Hamiltonian, i.e. 

= 7n7« where Hf = E„ e «7i7n- As long as only 
the zero-energy states contribute to $, the GGE is, in 
this case, spanned by the pair of MF. Thus the order 
parameter calculated within the GGE takes an elegant 
analytic expression, being given by the square modulus 
of the overlap between the (real) MF of the initial Hamil- 
tonian and those of the final Hamiltonian 



GGE 



E 



C/2 



dxu^{x)u { J\x) 



(10) 



where i and / denote the MF wavefunctions calculated 
with respect to Hi and Hf respectively. As shown in Fig. 
[T]the long-time average of $ after the sudden quench is 
in very good agreement with the corresponding quantity 
calculated within the GGE. For the TT -4 T quench the 
GGE predicts 'I'gge = since Hi has no MF; accordingly 
the real-time simulations provide values of §{t) close to 
0, see left panel Fig. [T] In the case of the T— >T quench 
the agreement is even more remarkable, see dashed line 
in right panel. Indeed the oscillations of <fr(t) (sizable 
only for Af < 100) are exactly centered around the value 
^gge- Thus we can infer that the GGE provides accurate 
predictions of the asymptotic quench dynamics of the 
Majorana order parameter. 

In Fig. [5] we illustrate the robustness of the T phase 
after a T — > T quench. The value $gge is plotted as 
a function of the strength of the initial magnetc field 
Vz l \ for a fixed value of the target V^K We see that if 
A < Vz^ < Vz^ the T phase is readily corrupted, even 
for small differences between the strengths of the initial 
and final magnetic fields. This happens in both the con- 



4 



1.0 
0.8 
0.6 

D 

0.4 
0.2 
0J 





N=50 










j-=10 2 


u 




— T=10 3 






— T=10" 






— T=10 5 






_ 





.0 0.2 0.4 0.6 0.8 1.0 
t/T 

0.020 , , , , , , , 

7V=200 



gO.010 



0.000 

1.0 0.4 0.5 0.6 0.7 0.8 0.9 1.0 




t/T 



FIG. 3: $ and lowest energy E op t as a function of time 
for TT — > T transitions wit different ramp duration and wire 
length. The protocol parameters are VP = 0.01 and vj f) = 
0.3 and the rest of parameters are the same as in Fig. [1] 
The discontinuities in &(t) observed for r = 10 2 are due to 
dynamical level crossings in the set {e w (£)}[2lj|. 



tinuum and in the lattice models. HVz > V z , instead, 
the continuum model predicts a significantly more stable 
T phase, since the system can afford larger changes in V z 
still maintaining values of $ close to 1 (see Fig. [5] dashed 
lines). 

Ramp protocols. — We now investigate the possibility 
of generating/maintaining the T phase after an arbitrary 
change in the strength of the magnetic field by consid- 
ering ramp protocols of finite duration r. In Fig. [3] we 
plot $(t) for different TT — > T transitions by varying r 
and J\f. It is seen that for sufficiently slow ramps the 
topological phase transition takes place and MF fully de- 
velop, in agreement with the adiabatic hypothesis. We 
mention, however, that the validity of the adiabatic the- 
orem is not obvious in the present contest, since the final 
Hamiltonian is such that its square Hj has a doubly de- 
generate ground-state, i.e. the pair of MF[27]. In the 
lower-right panel we also plot E opt (t). For J\f — 200 and 
t = 10 5 the formation of the MF is clearly visible as E opt 
approaches zero when V z (t) overcomes A (i.e. when Q 
becomes negative). Still, our results cleary indicate that 
the duration r required to reach $ = 1 increases with in- 
creasing size. For instance r ~ 10 3 is required to create 
MF in a wire of length J\f = 50, whereas r increases of 
about two orders of magnitude by enlarging the system 
to JV = 200. Interestingly this trend is reversed if one 
considers T — > T transitions. In this case the larger is 
the size of the wire, the shorter is the ramp-time needed 
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FIG. 4: Ramp-duration f required for a T — > T transition 
to reach the asymptotic value $ = 0.95 (boxes) and $ = 
0.9 (triangles) as a function of the length N ■ The protocol 
parameters are V z = 0.2 and V z = 0.3, and the rest of 
parameters are the same as in Fig. [T] 



to maintain the Majorana character after the transfor- 
mation. In Fig. [4] we show the duration f required for 
T — >■ T transition to reach a desired value $ of the or- 
der parameter (averaged over long times) close to 1. It 
appears that f is an overall decreasing function of TV 
(although with relatively small oscillations) that for the 
chosen parameters saturates to the finite value f w 80 for 
$ = 0.95 and f « 60 for $ = 0.9. We have checked (not 
shown) that this qualitative behavior does not change by 
varying the features of the protocol, thus providing an 
explicit measure of adiabaticity for protocols aiming at 
preserving MF. 

Conclusions. — We have studied the temporal evolu- 
tion of MF in driven ID quantum wires after the varia- 
tion of the system parameters according to different pro- 
tocols. In the case of sudden quenches, thermalization 
breakdown is observed, and the long-time behavior of the 
Majorana order parameter $ is well described within the 
GGE. Remarkably we find that the relaxed state does not 
display fully developed MF, no matter the initial Hamil- 
tonian is: for a TT— s-T quench the MF are not generated 
and $ remains close to zero, while for a T — > T quench 
the MF initially present get readily corrupted as the sys- 
tem is driven out of equilibrium. The adiabatic theorem 
is, instead, recovered for extremely slow ramp protocols. 
In the case of a slow TT —J> T ramp the system under- 
goes dynamically the topological phase transition within 
a timescale that increases by increasing the length of the 
wire. Finally we have provided an explicit estimate of 
adiabaticity to preserve MF during T — > T protocols, by 
pointing out the existence of precise temporal constraints 
relevant for the topological quantum computation. 
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